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ABSTRACT 


The results of seven different tests of a single stage 
axial turbine with converging-diverging stator nozzle are 
reported. From measurements in a special test rig, losses 
occurring in the stator and rotor blade rows were separately 
ealculated and the performance of the stage was also deter- 
mined. The rotor speeds varied from 9,500 to 18,600 r.p.m. 
and the pressure ratios varied from 1.75 to 3.25. The Mach 
numbers at the stator exit varied from 0.79 to 1.38. The 
results for this turbine are appraised and a procedure is 
demonstrated for smoothing loss coefficient data from the 
turbine rig. Test rig improvements reported anes the 
design and construction of a new flow nozzle and the revision 
of the data reduction programs to access a Hewlett-Packard 


Model 9867B Mass Memory unit. 
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ire eRODUCr EON 


The transonic turbine test rig in the Turbopropulsion 
Laboratory at the Naval Postgraduate School was designed 
to determine the effect of blading design on turbine effi- 
ciency and to allow the separate determination of rotor and 
stator losses. 

Preceding the work of Solms [Ref. 5] in which several 
different turbine geometries were tested, the separation of 
the losses in the blade rows was not achieved reliably. 
Progressive improvements in the hardware and instrumentation 
of the turbine test rig following each series of tests, and 
finally the simplification of the data reduction program, 
resulted in the successful separation of losses for particular 
configurations reported in Reference 5. However, anomalies 
remained to be explained, and the values obtained for the 
stator and rotor loss coefficients were scattered. 

The questions raised by the preceding test program mainly 
involved the results of the turbine configuration designated 
in Reference 5 as Turbine C. Turbine C had converging- 
diverging stator passages in an axial entry, single impulse 
stage which was designed to operate at a pressure ratio of 
3.71. Results obtained for this geometry discharging to 
atmosphere did not agree with results obtained at reduced 
pressure levels when plotted in terms of similarity varia- 


bles ("referred" quantities). Furthermore, the values 
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obtained for flow rate were of questionable accuracy since 
the existing flow nozzle was too large. Also the correction 
for a labyrinth leak rate became a more significant fraction 
of the turbine flow rate. 

The goal of the present work was to determine the blade 
row and stage performances of Turbine C, while resolving the 
anomalies reported in Reference 5. This would allow the use 
of the test rig to investigate effects of geometrical changes 
(stator-rotor separation and tip clearance for example) on 
the losses. The first step was to design and construct a new 
flow nozzle for the range of flow rates required by Turbine 
C. The flow nozzle was carefully calibrated against three 
different nozzles operating choked. The leakage rate for 
the labyrinth seal was also measured and a semi-analytic 
representation of the leakage rate as a function of pressure 
ratio was obtained which should apply at different tempera- 
tures. These developments are reported in Appendix A. The 
existing data reduction program [Ref. 5] was revised to 
access a Hewlett-Packard Model 9867B Mass Memory unit. The 
revised program is described fully in Appendix B. 

In the results of the test program given here in Section 
III, the effect of pressure level reported in Reference 5 
was not observed. The scatter in the loss coefficient however 
was not removed by the improvements in the accuracy of the 
measurements. In the discussion in Section IV, a method 


for smoothing the loss coefficient data is demonstrated 


uz 





which will then allow a useful study of the effect of axial 


and tip clearances to be carried out. 
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Pe elon UNG toot RIG INSTALLATION 


ee DESCRIPTION 

The test installation consists of three major components: 
an Allis Chalmers twelve stage axial flow compressor, an 
exhauster assembly, and the turbine test rig (TTR) itself. 

The compressor is the source of driving air for the TTR 
and for the exhauster assembly. Fig. A-1l shows the piping 
arrangement. Turbine air passes through the first settling 
tank into an eight-inch pipe containing a flow nozzle, into 
the second settling tank and into the turbine. 

Fig. 1 shows the plenum, the floating stator assembly, 
the rotor, and the dynamometer [Ref. 5]. Pressure ‘ratios of 
6:1 can be achieved when the system is hooded. The hood 
was needed to achieve near deSign pressure ratios in the 
tests reported here. Fig. 2 shows the turbine blading 
arrangement of the stator and rotor. Ref. 8 contains detailed 
description of the test rig hardware. 

The floating stator assembly depicted in Fig. 3 permits 
measurements of the stator torque while axial and rotational 
movements are detected by calibrated force transducers that 
are heat insensitive. These measurements allow the deter- 
mination of the axial and tangential velocity components 
at the stator exit [Appendix A, Ref. 5]. 

In this report one configuration designated Turbine C 


was tested, the geometry of which is shown in Fig. 4; Table 
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1 describes the geometry quantitatively. The stator blade 
profile is shown in Fig. 5. The blades generate a converging- 
diverging nozzle shape. Pressure measurements are taken at 
the locations shown in Fig. 5. The pressures necessary to 

the analysis of stator forces (Appendix B) are taken at the 


locations shown in Fig. 6 


B. TEST MEASUREMENTS AND ACCURACY 

MASS FLOW RATES: Appendix A gives a detailed analysis 
of the turbine flow rate and the labyrinth seal leak rate, 
and the results of calibration measurements. The results 
for the discharge coefficient are given in Fig. A-6. Fig. 
A-10 shows the Kinetic Energy Factor (Appendix A) obtained 
in labyrinth seal leak rate calibration and its eae neeee 
polynomial. The equations used in calculating the flow rate 
from measurements are summarized in Table A-1 and Table A-2. 

FORCES, TORQUES, TEMPERATURES AND PRESSURE: References 
8 and 9 give calibration procedures for the TTR. Identical 
procedures were employed here. Table II of Ref. 5 gives the 


expected accuracies of these measurements. 


C. TESTING AND DATA REDUCTION 

The TTR data collection system is described in Ref. ll. 
Appendix D in Ref. 5 gives a detailed explanation of the 
turbine test procedures. The data reduction program developed 
in Ref. 5 was revised to use the mass memory system. The 


revised program is described in Appendix B. 
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Seven tests were conducted of Turbine C at pressure 
ratios from 1.75 to 3.25. The parameters and conditions 


of the tests are summarized in Table II. 
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ae =LNTRODUCTION 

The reduced data from seven tests of Turbine C are 
listed in Table V to XI. Table IV gives a listing of the 
symbols and meanings of the Table headings used in Table V 
to XI. Plots of the reduced data are presented in Figures 
7 to 16. Table III gives an explanation of symbols used in 
these figures. 

Tests 1 and 2 were conducted with the hood off. Tests 
3, 4, 5, 6 and 7 were made with the hood installed. Tests 
5 and 6 were repeats of tests 3 and 4 respectively. The 


results are discussed in the following paragraphs. ~ 


eee TORQUE 

Fig. 7a and 7b give the results of the referred torque 
versus referred speeds at a constant pressure ratio. A 
nearly linear dependency exists for all seven runs with 
consistency in performance for the 1.75 and 2.25 hooded and 
unhooded pressure ratios. 

The results from the unhooded tests (Fig. 7a) agree with 
the results of the hooded tests (Fig. 7b), at the same 


pressure ratio, to within 1 percent. 


See BFFICIENCY (TOTAL-TO-STATIC) 
Figures 8a and 8b give the results of the total-to-static 
efficiency versus referred R.P.M. Figures 9a and 9b show 


plots of efficiency versus isentropic head. The maximum 


ey, 





efficiency occurred at similar values of the isentropic 
head for all pressure ratios. Ref. 17 relates efficiency 
to isentropic head. The peak efficiency was achieved at an 
isentropic head coefficient between 4.0 and 4.5 for both the 
hooded and unhooded tests. 

The efficiency is seen to increase as the pressure 
ratio is increased, approaching 80 percent at the pressure 
ratio of 3.25. The efficiency would be expected to increase 
to the design pressure ratio of 3.71, at which the stator 


nozzles should be just correctly expanded. 


D. THEORETICAL DEGREE OF REACTION 

Figures 10a and 10b give the results of theoretical 
degree of reaction versus isentropic head coefficient at 
different pressure ratios. Ref. 5 and Ref. 17 give the 
equation for the theoretical degree of reaction. For all 
runs the theoretical degree of reaction was negative at 


lower speeds. For a negative value, P./P must be greater 


=O 
than Pi/PLo: The values become poSitive only at high r.p.m. 
in both the hooded and unhooded configurations, passing 
through the design value of zero at a value of the isentropic 
head coefficient corresponding to the maximum in efficiency. 
Thus, negative values of the theoretical degree of reaction 
occur at speeds lower than the design speed at a given 
pressure ratio. Typically, there is seen to be also a loss 
of three to five percent in efficiency from the maximum 


value as the isentropic head coefficient is increased 100 


percent above the design value. 
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Results from hooded and unhooded tests are seen to 
agree with the exception of Run 3 in Figure 10b for which 
the values are low. Lack of agreement for the data of this 
test was found in all but the overall stage performance 


measurements. 


Bm. EFFECTIVE DEGREE OF REACTION 

Figures lla and llb give the results for the effective 
degree of reaction versus isentropic head coefficient. At 
all pressure ratios, hooded and unhooded, the effective 
degree of reaction was negative at lower speeds. Ref. 17 
states that the effective degree of reaction is a means for 
judging whether the flow is being accelerated or decelerated 
in the rotating cascade. The negative degree of reaction 
in Fig. lla and 11b show that the flow is decelerated at 
lower r.p.m. and becomes accelerated as r.p.m. is increased 
(Isentropic head is decreased). In no case should the design 
value of the effective degree of reaction be less than zero. 
As stated above, operation at speeds less than design speed 
results ina loss in stage efficiency, partially due to the 
undesirable deceleration in the rotor. 

The results for Run 4 (at a pressure ratio of 3.25) in 
Figure llb reflect the marked improvement in efficiency at 


pressure ratios approaching design. 


iE. STATOR LOSS COEFFICIENT 
The results for the stator loss coefficient are shown 


in Fig. 12a and Fig. 12b. While there is some scatter in 


ie 





the data, the results from hooded and unhooded tests are 
seen to be in reasonable agreement. There is a slight 
increase in the stator loss as speed 1s increased at any 
pressure ratio. The stator loss for the highest pressure 
ratio (3.25) is significantly lower (13-17 percent) than for 
the lower pressure ratios (18-25 percent). This reflects 
the poor performance of a converging-diverging nozzle at 


pressure ratios less than design. 


Gene ROTOR LOSS COEFFICIENT 

The results for rotor loss coefficient are shown in 
Figure 13a and Figure 13b. For completeness, the results of 
Run 3 have been included where the rotor loss coefficients 
were evaluated to be positive. Negative loss coefficients 
were obtained at higher speeds in Run 3, indicating an 
error in at least one element of the primary data. 

The rotor loss coefficient is the performance parameter 
most sensitive to small errors (and possibly unsteadiness) 
in the measurements. The sensitivity arises because it 
involves most of the parameters (see Ref. 5) previously 
calculated in the form of a ratio of small differences. In 
view of this, if the results of Run 3 in Fig. 13b are excluded 
from consideration, the rotor loss coefficient for the tests 
are reasonably consistent. The results for hooded and 
unhooded tests agree to within 20 percent at corresponding 
pressure ratios. There is also a definite decrease in the 


rotor loss coefficient as pressure ratio is increased toward 
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the design value, and a consistent decrease as the speed 
is increased at a fixed pressure ratio. The simultaneous 
reduction observed in the rotor and stator loss coefficients 
at higher pressure ratios is consistent with the increase 


measured in the stage efficiency (Fig. 8 and Fig. 9). 


H. FLOW ANGLES 

An example of the variation in the flow angles at the 
stator and rotor exits as a function of the isentropic 
Mead coefficient is shown in Fig. 14. As was found in all 
cases tested in Ref. 5, the variations in stator exit flow 
angle and rotor exit relative flow angle are very small. 
Also, the total variation of the relative flow angle into 
the rotor is only about 6 degrees for a 50 percent change 
in r.p.m. The results in Fig. 14 are qualitatively repre- 
sentative of the results for each of the seven tests. The 
exact variations in the measured angles can be seen in 
Tables V to XI. 

The stator exit flow angle is shown in Fig. 15a and 
Fig. 15b as a function of isentropic head coefficient for 
the seven tests. It can be seen that the measurements are 
in the range of 75 to 77 degrees, with the lower value 
occurring at the highest pressure ratio (3.25). It can be 
seen in Fig. 5 that the suction surface of the stator blade 
is designed to exit at 75 degrees, and the pressure surface 
at a smaller angle. Thus the measured stator exit angle 


mmedrs tO be Slightly larger than would be expected; however, 


zee 





the measured values when the nozzle passages are nearly 
correctly expanded approach the geometrical angle. 

The large values of stator exit flow angle at lower 
pressure ratios, and the higher values of stator loss 
coefficient are probably the results of the presence of flow 
separation. The effect of flow separation would be reflected 
in two ways. Firstly, the flow angle of the separated flow 
could depart significantly from the trailing-edge angle. 
Secondly, the existence of a region of separated flow could 
seriously affect the averaging process which is implicit 


in the data reduction. 


ie FLOW RATE 

The referred flow rate is shown as a function of the 
isentropic head coefficient for the seven tests in Fig. 
l6a and Fig. 16b. The referred flow rate was constant at 
1.02 lbs/sec to within the accuracy of the measurements 
for all speeds and all pressure ratios when operating 
hooded (Fig. 16b). During unhooded operation (Fig. l6a), 


the results gave a constant value of 1.01 lbs/sec. 
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A. RESOLUTION OF PREVIOUS ANOMALIES 

In the results presented above, the difference in the 
turbine performance operating hooded compared with unhooded, 
1s seen to be small. This is in contrast to the results pre- 
sented in Reference 5. Three factors contributed to the 
resolution of this anomaly; the accurate measurement of the 
turbine flow rate, installation of temperature insensitive 
load cells within the hood to measure stator loads, and 
corrections that were made to the data reduction program. 
There remains to be explained a difference of 1% in the 
flow rate between hooded and unhooded operation. More data 
are needed from unhooded tests to confirm that this differ- 
ence is repeatable. Because of the sensitivity of the perfor- 
mance parameters to the flow rate, the residual disagreement 
between hooded and unhooded results might be resolved if 


the difference in flow rates were explained. 


B. FLOW. RATE 

The use of a choked nozzle to calibrate a large flow 
nozzle was successful (Appendix A). With the blockage 
factor correction applied, the results for different sizes 
of the choked nozzle agreed in the range for which the flow 
haces overlapped. 

The labyrinth leak rate calibration measurements were 


well represented by the analysis given in Reference 7. The 
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value obtained for the Kinetic Energy factor was higher than 
values given in Reference 5; however, the geometry of the 
turbine testing labyrinth was not known exactly. There was 
some eccentricity in the alignment of the stator section 
within the gland, so that the clearance was nonuniform to 
an unknown degree. 

The flow rate through the turbine, when expressed as 
a referred quantity, was measured to be a constant over the 
range of pressure ratios tested. (However, there was scatter 
in the measurements that resulted from unsteadiness in the 
measurements of nozzle pressure drop.) Since the stator 
was choked, the constant value was expected. The value 
of the constant was consistent with the flow rate computed 
for the stator using a calculated value for the stator 
throat area. A value of the blockage factor within 1% of 
unity was required to be applied to the geometrical stator 
area to obtain total agreement. Since the exact area of the 
throat where choking occurs in the stator is not known to 
this accuracy, it is concluded that flow nozzle and stator 


throat measurements of the flow rate are in agreement. 


fee LOSS COEFFICIENTS 

The scatter observed in the results for the loss coeffi- 
cients makes a comparison between predicted and measured 
values inconclusive. There are trends in the data which are 
discernible, the rotor losses are consistently higher than 


the stator losses, for example. However, the scatter does 
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not allow the effect of parameter changes on the separate 
loss coefficients to be determined. The scatter is the 
result of the extreme sensitivity of the loss coefficient to 
the separate measurements on which they depend. For 


example, the stator loss coefficient Ge) 1S given by 


Z 


= = PEM pe aces 
Cs ae y-1 (2p) 
= Y 
ip Py 


where Xy is the non-dimensional velocity and Py 1s the 
calculated pressure ratio (p,/P,.) at station 1. Differ- 


entiation of Eq. (1) gives, if Xy 1s constant, 


y-L 
aie ce 2 Deny dP 
ao -( (8) (+) (+ _t (On 
Cs Gs 1 Y )]J an 
1 Bh Y 


Baking the sonic condition of P, = 0.532, with y = 1.4 and 


a 
= = 0.2, the quantity in brackets has the value 5.78. Thus 
a 1% error in calculating Py produces -6% error in the stator 
loss coefficient. Since Py is derived (as is X,) from 
measurements scanned over a period of more than a minute, 
with some variation in operating point (as well as an 
observable fluctuation in the flow nozzle pressure drop) 


occurring during the data scan, the scatter in the loss 


coefficient data is understandable. 
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Before the loss coefficient can be used as a measure 
of performance, a method of smoothing the data while 
retaining accuracy is needed. Two attempts were made 
here. First, an attempt was made to smooth the calculated 
values of the interstage pressure (Dy), by comparing Py 
to the pressures measured at the hub and tip. Figure 17 
Shows the variation Py and the measured pressure with 
isentropic head coefficient for one particular test. 
Clearly, the calculated average pressure at station 1 must 
be smoothly behaved if the hub and the tip pressures are 
smoothly behaved. An attempt to derive a smoothing function 
for the interstage pressure is shown in Fig. 18, where the 
data from several tests at the same pressure ee nealley « elek= 
plotted in a dimensionless form. The figure clearly shows 


the degree of scatter in the values of P On the basis of 


1° 
these data above, however, the use of a Single polynomial 


mer the quantity (Py-PL p/P ) as a function of 


tip hub 
isentropic head coefficient cannot be justified. When a 
polynomial was derived from the data of a Single test, and 

the loss coefficients for the stator and rotor were recomputed, 
the results shown in Fig. (19a) and (19b) were obtained. 

The second attempt to smooth the loss coefficient 

data involved effectively the elimination of scatter in the 
input data from which the interstage pressure and velocity 
were calculated. The referred flow rate was determined to 


be constant when the stator was choked. By setting the 


referred £low rate equal to the constant value (1.02 or 1.01, 
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hooded and unhooded, respectively) and continuing the data 
reduction as before, the results shown in Fig. 20a and 20b 
for the loss coefficient were obtained. A considerable 
reduction in scatter 1S evident, with the major benefit in 
the rotor loss coefficient. Further data must be examined 
to fully appraise the benefit of this smoothing technique. 
A third technique, so far not attempted, is to calcu- 
late Py in the data reduction before calculating the axial 
velocity component at station 1. The values of Pi could 
be smoothed using a polynomial curve fit and the smoothed 
values of Py used to calculate the velocity. An application 
of all three techniques should produce loss coefficient 


distributions which are smooth and yet accurate. 


aa 





V. CONCLUSIONS AND RECOMMENDATIONS 


The following conclusions are made for turbine configura- 
mon C. 

1. The loss coefficients are extremely sensitive to 
input data, especially the rotor loss coefficient. Repeata- 
bility in calculating the loss coefficient directly from 
measurements was poor. 

2. Methods were demonstrated to smooth the loss coef- 
ficient data without sacrificing, and possibly improving, 
the accuracy. 

3. For all pressure ratios there was a velocity 
deceleration through the rotor and a pressure rise ‘across 
the rotor at lower r.p.m., indicating perhaps a shock wave 
mewoOr in the rotor. 

The following recommendations are made: 

1. Apply the three data smoothing techniques described 
above to all data from Turbine C in sequence: input constant 
referred flow rate, calculate interstage pressure and smooth 
using measured hub and tip pressures, then calculate 
interstage velocity. 


2. Determine a more accurate input for K the ratio 


5 
of the momentum average velocity to the mass-average velocity 
and examine the effect on the calculated stator exit flow 


angle. 
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3. Program the quadratic curves for the theoretical 
loss coefficient in references 16 and 18 in order to reduce 
the input error. 

Finally, it is concluded that following the recommenda- 
tions listed here, the effect of parameter variations on 
the blade row losses could be examined satisfactorily using 


the turbine test rig. 
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TABLE IV 


DEFINITION OF COLUMN HEADINGS IN TABLES V_TO XI 


SYMBOL 


ve 


V2 


VAL 


VA2 


Al 


A2 


BI 


B2 


ZS 


ZSTH 


ZR 


2RTH 


ZR* 


ZL 


DEFINITION 
VELOCITY AT STATOR EXIT PLANE (FT/SEC) 

VELOCITY AT ROTOR EXIT PLANE (FT/SEC) 

AXIAL VELOCITY AT STATOR EXIT PLANE (FT/SEC) 
AXIAL VELOCITY AT ROTOR EXIT PLANE (FT/SEC) 
TANGENTIAL VELOCITY AT STATOR EXIT PLANE (FT/SEC) 
TANGENTIAL VELOCITY AT ROTOR EXIT PLANE (FT/SEC) 
MACH NUMBER AT STATOR EXIT PLANE 

AXIAL MACH NUMBER AT STATOR EXIT PLANE 

MACH NUMBER AT ROTOR EXIT PLANE 

AXIAL MACH NUMBER AT ROTOR EXIT PLANE 

FLOW ANGLE AT STATOR EXIT PLANE (DEGREES) 

FLOW ANGLE AT THE ROTOR EXIT PLANE (DEGREES) 
RELATIVE FLOW ANGLE AT ROTOR INLET PLANE (DEGREES) 
RELATIVE FLOW ANGLE AT ROTOR EXIT PLANE (DEGREES) 
STATOR LOSS COEFFICIENT 

THEORETICAL STATOR LOSS COEFFICIENT 

ROTOR LOSS COEFFICIENT 

THEORETICAL ROTOR LOSS COEFFICIENT 

ROTOR CARRY OVER LOSS COEFFICIENT 

ROTOR INCIDENCE LOSS COEFFICIENT 

PRESSURE RATIO Po! Po 

ISENTROPIC HEAD COEFFICIENT 

REFERRED RPM (R.P.M.) 


2 





WREF 


ETA 


RTH 


REFF 


Pl 


EZ 


PS 


P4 


p> 


P6 


By 


P8 


Ee? 


KBLOCK 


STPR 


RPM 


REFERRED FLOW RATE ( LBM/SEC ) 

TOTAL TO STATIC EFFICIENCY 

REFERRED HORSEPOWER ( H.P. ) 

THEORETICAL DEGREE OF REACTION 

EFFECTIVE DEGREE OF REACTION 

PRESSURE RATIOS CORRESPONDING TO TAP LOCATIONS 
INSEEGURE S- ERESSURES ARE REFERRED TO STATOR 


INLET TOTAL PRESSURE 


STATOR BOCKKAGE FACTOR 
PRESSURE RATIO ACROSS THE STATOR P /P, 
O 


ROTOR SPEED CR, P.M.) 
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APPENDIX A 


FLOW RATE DETERMINATION 


A-l INTRODUCTION 

A precise measurement of the flow rate is necessary 
if the performance of an individual blade row is to be 
determined accurately from measurements made with the turbine 
test rig. Anomalies that occurred in the results from the 
first tests of the present Supersonic turbine (Turbine C 
in Ref. 5) were possibly due to inaccuracies in the speci- 
fied flow rates. Since the flow rate for Turbine C was 
smaller than for the turbines previously tested in the rig, 
small pressure differentials were measured with the existing 
nozzle. Also, the higher pressure ratios at which Turbine 
C was operated resulted in increased leakage through the 
labyrinth seals, the geometry of which might have change 
slightly on reassembly. Therefore, before new tests of 
Turbine C were begun, a new flow nozzle was provided and 
careful measurements were made of the labyrinth leak rate. 
The flow nozzle design and calibration are described in 
Section A-2. The labyrinth leak rate determination is given 


in Section A-3. 


A-2 MEASUREMENTS OF THE TOTAL FLOW RATE 
A-2.1 FLOW NOZZLE DESIGN 
The flow nozzle is positioned in a pipe through which 


air is supplied to the turbine (Fig. A-1l). The nozzle 
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generates a differential pressure to which the flow rate 
is related. 
For a one-dimensional steady flow, the mass flow rate, 


W, iS given by 
W =o AV (A=1) 
where 9 is the density, V is the velocity, and A is the 


cross sectional area. 


Using the perfect gas equation of state 


0 = — (A-2) 
and Bernoulli equation 
P.-P=59Vv | (A-3) 
Eq. (A-1) for low velocities becomes 
: 2P, 
; W=A Rr, (Pt - P) (A-4) 


The symbols P and T are pressure and temperature, respec- 
tively, and the subscript t denotes stagnation conditions. 
This equation, with explicit corrections for compressibility 
and thermal expansion, is used to calculate the flow rate 


through a flow nozzle from measurements of temperature and 


To 





pressures. Withe the addition of the correction factors and 
with conversion factors added to account for the units of 


measurements [Ref. 1], Eq. (A-4) becomes 








Wy = pon les8s Dee Ay Ky ¥, [—G— (1bs/sec.) —(A-5) 
¥2.036 te 
where 
Du = Diameter of nozzle (inches) 


Ay = Thermal expansion coefficient 
Y, = Compressibility coefficient 
P,, = Pressure at the nozzle (psia) 


Ah = Water differential (in. H,0) 


T, = Temperature (°R) 
0.16384 _ Conversion factor 


72.036 


and lastly, K,, is the discharge coefficient. 


N 
The discharge coefficient accounts for the fact that the 
flow is not one-dimensional, and that the pressures measured 
are not Py and (Py - P). Generally, the discharge coefficient 
has a value above 0.95 for a standard ASME nozzle at high 
Reynolds numbers. Establishing the relationship between 


Reynolds number and the discharge coefficient is what con- 


stitutes the calibration. 


80 





The flow rates expected in the supersonic turbines were 
1.3 to 3.5 lbm/sec. A flow nozzle with a diameter of 3.25 
inches was chosen (by applying Eq. (A-5) with unity for the 
coefficients) to give sutiable water column measurements of 
approximately 9.8 to 58.1 inches. The nozzle design followed 
the criteria given in Ref. 1 for "Low Beta series" designs. 
The geometry of the flow nozzle is shown in Fig. A-2 and 
views are given in Fig. A-3. The locations of the pressure 
taps are shown in Fig. A-2. Four throat taps are located 
1.5 nozzle diameters from the nozzle face and are manifolded 
together. The upstream pressure 1s measured uSing taps in 


the upstream flange. 


A-2.2 FLOW NOZZLE CALIBRATION 
The discharge coefficient is given in terms of measure- 
ments, if the flow rate is determined by independent means. 
One of the most accurate methods of determining the flow 
rate of a gas is by choking the flow at a known area. This 
1s because the mass flux is nearly uniform where choking 
occurs, and because the boundary layer can be kept thin in 
a rapid smooth contraction from a relatively large pipe. 
Hence, in order to determine the discharge coefficient of the 
new flow nozzle as a function of Reynolds number, the 
arrangement shown in Fig. A-4 was used. 
A-2.2.1 Method 
Flow nozzle calibration runs were made at con- 


trolled supply pressures from 15 psia to 45 psia. Three 
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different orifice plates with diameters of 1.6, 2.065, and 
2.24 inches were mounted in turn to provide choked flows 
over the desired range of Reynolds numbers. The differen- 
tial pressure at the flow nozzle was varied in increments 
of 2 to 4 inches of water. The range of flow rates over 
which the orifice plates were choked at the supply pressure 
was from 1.33 to 3.55 lbm/sec. The maximum pressure ratio 
obtained was 2.92. 

A-2.2.2 Analysis 

The discharge coefficient was given by Eq. (A-5), 
and the Reynolds number was determined by 

48 Wy 8 


Rey = . (A-6) 


t D 
p u 


where u is the viscosity and 8 is the ratio of the diameter 
of the nozzle to the diameter of the pipe. The flow rate, 
which was the same for the nozzle as for the choked plate, 


was determined by the following equation 


where 


tg 
I 


pressure at choked nozzle (psia) 


area of choked nozzle (in?) 


yp 
il 
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~~ 
I 


blockage factor 


mW 
il 


gas constant (ft-lb,-/lbm - °R) 


re 
il 


Temperature at choked nozzle (°R) 


The function [ is a known function of the specific heats and 


is given by 


z moe eee 2 


The blockage factor K, accounts for the boundary 


B 
layer and for the axisymmetric surface; Kp is given by 


fRef. 3] 
KK. =l-=> (A-9) 


where the star on D signifies sonic conditions, and §* 


N 
is the boundary layer displacement thickness. 
The boundary layer growth was estimated using 


the method of Falkner and Garner given in Ref. 6. The 


displacement thickness, 6*, is given by 


2 n 
* n+l —— 
fa Ha yrrt (A-10) 
Sj IL 
Re n+1 
1 
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where S) is the distance from the stagnation point to the 
choking station. H is the form factor, and a and n are 
empirical constants. wy is a factor which depends on the 
pressure gradient along the surface. Re, is the Reynolds 
number based on conditions at the throat and Sj- 

It is shown in Ref. 6 that the value of w for 


a Similar contraction was 0.38. With w = 0.38, and with 


H= 1.4, a = 0.0076 and n = 6 for air, Eq. (A-10) becomes 


s* = Oreos 4 z 
0.1429 “1 
Re 
i: 
The Reynolds number is given by 
Ww, - 
Re, = 12 — — (A-12) 
- ie 


where Wy 1s the viscosity at the throat where the tempera- 
ture is 5/6 times the stagnation temperature. 
For each operating condition for each orifice 


plate, Re, 


(A-11) 


was calculated from Eq. (A-12) using W calculated 


from Eq. (A-7) with K, = 1.0. 6* was then calculated using 


B 


Eq. (A-1ll) and K, obtained from Eq. (A-9). W., was recalcu- 


B N 


lated from Eq. (A-7). It was found that the blockage factor 


could be taken to be constant with pressure for each choked 


plate at the following values: 
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PLATE # THROAT DIAMETER Sy fans.) Kp 

as re aol 029928 
a ZA 065 SereLo OEI725 
ae 2.24 Bi 2 029904: 


The flow chart for the analysis of the calibra- 
tion test results is given in Fig. A-5 and the equations 
used are summarized in Table A-l. 

A=-2.2.3 Results 

The nozzle calibration extended over a range 
of Reynolds numbers from 0.081 x 10> O00. 2 Lb 30x ie. The 
results are given in Fig. A-6, where the data for three 
different orifices are shown for the range of presstres 
for which the flow was choked. A second order polynomial 


was found to represent the flow nozzle coefficient as a 


function of the Reynolds number. 


Rey Rey 2 
Ky = 1.0272 - 0.1598 cama? + 0.3805 era 
10 10 


(A-13) 


where Rey 1s given by Eq. (A-6), was used to reduce data 
from tests of the turbine test rig. 

A-2.3 APPLICATION 

The discharge coefficient is a function of Reynolds 
number. However, in application the flow rate is to be 
determined and therefore the Reynolds number is unknown. 


The determination of the flow rate from measurements requires 


the solution of three equations which have the general forms: 
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W.. = a, K (A-14) 


(A= 1:5) 


x 
\ 
ms 

Z 
Cy 
0 

AL 


Re, = a, W a1) 


These equations represent Eq. (A-5), Eq. (A-13), and Eq. 
(A-6), respectively, where Oy and a. are known in terms of 
the geometry and measurements for any operating point. 
Clearly, by substituting for Rey and Ky MS ING sea) 
and Eq. (A-16) in Eq. (A-14), a single function for the 
unknown W 1s obtained. 

An iterative procedure is used to solve the three 


equations for the unknown flow rate. The procedure is 


given in Appendix B. 


A-3 MEASUREMENT OF THE LABYRINTH LEAK RATE 

Reo sbabyrinth LeakeRate Calibration Test 

Since the leak rate of the plenum labyrinth cannot be 
measured during turbine operations, the relationship of the 
leakage flow rate to supply and hood pressures and supply 
temperature was determined in a separate calibration test. 
For this test, the main flow through the turbine was blocked 
by a plate and the labyrinth was allowed to leak into the 
hood as in normal operation. The flow was passed through a 
two-inch pipe (Fig. A-4b & Fig. A-1), containing an ASME 


standard sharp-edged orifice. The usual supply pipe was 
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blocked. The flow was discharged through the labyrinth 
seal into the hood in which the pressure was controlled 
by the ejector. 

In determining the leakage rate using the ASME orifice, 
vena contracta taps were used. The differential pressure 
was measured in inches of water on a manometer which had 
O.l-inch graduations. The upstream pressures in the lab- 
yrinth and inthe nozzle were measured in inches of mercury 
on U-tube manometers with atmospheric pressure as a refer- 
ence. The supply pressure was varied from 14 to 83 psia, 
while the differential pressure across the vena contracta 
taps varied from 2.75 to 37.3 inches of water. The maximum 
pressure ratio across the labyrinth was 5.6. 

A-3.2 Calculation of the Leakage Flow Rate 

Since the flow orifice was in a standard installation 
according to the criteria given in Ref. 1, Table 5 of Ref. l 
was used to obtain the discharge coefficient as a function 
of Reynolds number. In order to program the data reduction, 
the data of Table 5, for the appropriate values of Beta, 
were first approximated in the range of Reynolds numbers 
from 8,000 to 1,000,000 by a polynomial as shown in Fig. 
A-7. The discharge coefficient at each operating point was 


then given by 
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5 
O2o02098259°+ 0,004688755 (=o) CAS) 


K = 
N N 
-310°.2 oe 
~ 1.250025 x 10° (32-)* + 2.0320771882 x 1077 (22°) 
Re Re 
N N 
-5 10°.4 -7 10°.5 
- 1.56980769 x 10 (——) + 4.5277098 x 10 ==) 
Rey Rey, 


The flow rate was calculated at each point using Eq. 
(A-5), using the expressions for the compressibility and 
thermal expansion coefficients given in Table A-2, and the 
discharge coefficient given by Eq. (A-17). 

A=3.3 Analysis of the Results 

A-3.3.1 Simple Method 
A non-dimensional flow rate, or flow rere 


for the flow through the labyrinth can be defined as 


W, ¥ (R/g) Tr. 
6 = (A-18) 


ei 

During the calibration tests, the stagnation 
pressure (P,) wasS varied in steps for constant values of the 
hood pressure (P)). The stagnation temperature was approxi- 
mately constant. The flow function was calculated for each 
operating point using Eq. (A-18) and the values plotted as 
a function of the pressure ratio (P)/P,) - The results are 
given in Fig. A-8. | 

A second order polynomial was Found EO L£Ltethe 


data to an accuracy of + 4% as shown in Fig. A-8: 
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— -_ 2 = 
® = 0.2721 + 0.0365 {P,/P_) - 0.2357 (P,/P_) (A-19) 


While Eq. (A-19) represents the calibration tests 
results quite well, it is an empirical result for a limited 
Variation in test parameters. Operation at different 
temperatures, for example, could give a different behavior. 

Pao Uoemor scien kinettc ENnengy,) Factor 

It is shown in Ref. 7 that, under reasonable 
assumptions, the flow through a labyrinth with many teeth 


(N) can be represented by 
@ = 8 Ke (A-20) 


where 





2 
1-(P,/P,) 


i (A-21) 
N-Ln (P| /P,) 


and Ke is the "Kinetic Energy Factor". Here, the dependence 
of the leak rate on the upstream and downstream pressures 
(8) has been obtained analytically. The Kinetic Energy 
Factor is a measure of how much kinetic energy is recovered 
as stagnation pressure as the flow passes saccessSive teeth. 
It is expected to be a function of Reynolds number, and 
hence would include the correct behavior with stagnation 


temperature. 
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The geometry of the labyrinth is given in Fig. 
A-9. Using the notation in Fig. A-9, the Reynolds number 


for the leakage flow is given by 


12 We ec 
R = —-—— A-22 
mi AL Hy ( ) 
where 
AL = 2 1.) C (A-23) 


and Uy, is the viscosity based on the labyrinth supply 
temperature, Th: 
For the flow rates obtained in the calibration 
tests, 6 was obtained from Eq. (A-18). 8 was calculated 
uSing Eq. (A-21) and Ke determined using Eq. (A-20). The 
Reynolds number was obtained using Eq. (A-22) with Eq. (A-23). 
The Kinetic Energy Factor was plotted as a function of the 
Reynolds number and the result is shown in Fig. A-10. A 


fourth order polynomial was found to fit the data to an 


fecuracy of + 2%: 


Ke = 0.601373889 + 4.36695141 x 10% (5) 
10 
- 2.4530860 x 10°’ (—28)* + 6.5672233 x 10° 1+ (RS)? 
3 5 
10 10 
meAmacas 72. lon. hs (A~24) 
10 
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It is noted that the use of Eq. (A-24) and Eq. 
(A-21) gives a factor of two improvement in accuracy over 
the use of the purely empirical result given as Eq. (A-19). 
Also, the dependence of the Kinetic Energy Factor on the 
supply temperature is at least partially included in the 
representation as a function of Reynolds number. 

The flow chart for the calibration data reduction 
1S given in Fig. A-ll and a summary of the equations is 
given in Table A-2. 

A-3.2 APPLICATION 

In order to use Eq. (A-24) to calculate the leak rate 
at a given operating point, an iterative procedure is 
necessary Since the Reynolds number is again a function of 
the flow rate. The procedure 1S Similar to that followed 
in calculating the total flow rate (Section A-2.3) and is 


described in Appendix B. 


or 





TABLE A-l 


SUMMARY OF FORMULAS FOR NOZZLE CALIBRATION DATA REDUCTION PROGRAM 
0.16384 


2 
= Doe «oY 
[ 2.036 eee 1 ee 


EXPANSION 140.00193 
COEFFICIENT 


COMPRESSIBILITY LemOnOS 246 
COEFFICIENT 


x 0.06333 Ty 
VISCOSITY - 


TIS Ty * L 


FLOW 
FUNCTION: 


0.21364 


0.129 
Re, 
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TABLE A-2 


SUMMARY OF FORMULAS FOR LABYRINTH CALIBRATION REDUCTION PROGRAM 














; RI 
FLOW ¢ _ 2.036 Wy 3 
FUNCTION: 
FLOW p 0.16384 ae 
RATE: J 2.036 Te 
1 - Ch 
BETA: a citin owe 0 


REYNOLDS 





KINETIC $ 
ENERGY Ke 
FACTOR: 8 





THERMAL Ease 
EXPANSION ‘ 1+ 0.00193 | ————_— 
COEFFICIENT: : 100 


> 
I 


COMPRESSIBILITY 


COEFFICIENT: 








BETA: P = » Dy = 3.25 ins., D, = 7.975 ins. 


25 
1,153 x 10° x 0.063337 Ty 
VISCOSITY: UM = 





195 J 2a/ Ty soe 
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APPENDIX B 


TURBINE TEST RIG (TTR) DATA REDUCTION AND PROCESSING 


fe =LNIRODUCTION: 

This appendix documents the current data reduction 
program which was revised and updated from the program used 
in Ref. 5. The present program uses the "Mass Memory" disc 
storage unit attached to the Hewlett-Packard Model 9830A 
Calculator. The equations developed in Ref. 5 apply to 
the present program. 

A description of the programs is given in Section 
B-2 and the variables are listed and defined. Procedures 


for running the program are given in Section B-3. 


B-2 DESCRIPTION OF PROGRAMS: 

The reduction program is divided into eight separate 
sub-programs. Fig. B-l shows the contents of each sub- 
program. Table B-1l gives the channel and port numbers for 
each of the test measurements and the respective matrix 
elements assigned in the data reduction program. Table B-2 
defines the variables and equations (Ref. 5 & 16) used in 
the program. Table B-3 is a record or score sheet of the 


variables used in the program. 


B-3 PROCEDURE FOR DATA REDUCTION PROGRAM USING MASS MEMORY: 
ie) CHECK: 
Mass Memory to UNLOAD Position 


a.) 
b.) ROM installed 
c.) NO cassettes in transport 
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EO. 


. 


a2 . 


5 . 


A. 


HS < 


lULieonmmdr 9650A Calculator 

Turn on HP9866A Printer 

Turn on HP 9868A I/O Expander 

Turn on HP9863A Tape Reader 

Turn on HP 9867B Mass Memory 

Turn on HP 11305A Controller 

Load platter number TTR when the DOOR UNLOCKED light is 
ey 

Switch to LOAD position; wait 30 seconds until DRIVE 
READY light is lit. 

Key inse GET "TTR 1:"“ Wait for completion 

Press "RUN EXECUTE." 

The display will read "ENTER NEXT RECORD # ON DATA FILE." 
Input your next record number for raw data storage and 
press “EXECUTE." 

The display will read "TAPE: 1ST HOLE?-ON START?:-CONT." 
By pressing "CONTINUE, EXECUTE" the paper tape will be 
read. 

After the paper tape is read, the display will read 
"CORRECTIONS TO DATA? YES=1, NO=0;" This option permits 
you to correct any erroneous corrections. PRESS "1 
EXECUTE" if you have corrections. Press "0 EXECUTE" 

if no corrections; go to step 19. 

With Corrections to be made, the display will read 
"PRESS PRT ALL KEY FOR RECORD." This enables you to 


have a record of your data points. 
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1a 


a7 


eS. 


no 


20. 


21. 


The display will read "ENTER CORRECTION AS MATRIX 
ELEMENTS." When a scanivalve is detected to have given 
an erroneous value, the manometer board can be read 

and correct value recorded. See Table B-2 for the 


appropriate matrix identification. 


il 
J 
=e 


The display will read "ENTER CORRECT VALUE 
Input the correct value and press "EXECUTE." 

The display will read "ANY MORE CORRECTIONS? YES=1, 
No=0." Press "1 EXECUTE" if more corrections; go to 
step 17. 

1.5 inches of paper will emerge from the printer and 

the raw data will be printed. Afterwards the display 
will read "STORE DATE? ENTER YES=1, NO=0." If. raw 

data is to be stored in Mass Memory press "1 Execute." 
The raw data is stored. Below the tables will be 
printed the words "THE RAW DATA IS STORED IN RAW DAT 
FILE #? Along with the file number that was input in 

step 12. If data is not to be stored press "0Q EXECUTE" 
and the words "THIS DATA WAS NOT STORED" will be printed. 
The calculator will automatically CHAIN with "TTR 2." 
Temperatures of the nozzle, labyrinth and Tho will be 
printed in degrees Rankine. The labyrinth seal flow 
rate, its Reynolds number and kinetic energy factor 

will be printed. Also, the nozzle discharge coefficient, 


nozzle flow rate, its Reynolds number, and flow rate 


will be printed. 
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ZZ. 


23). 


24. 


BD. 


Zio). 


The display will read "CALC NEW W-DOT? YES=1, NO=0." 

If you elect to calculate a new theoretical flow rate, 
based on a stator blockage factor of 0.965, and employ 
this value in- the remainder of the program, press "l 
EXECUTE." However, if you use the original flow rate, 
press "0 EXECUTE." Go to step 24. 

The election for a new calculation of flow rate was 
chosen; a new flow rate, the stator flow function and 
blockage factor will be printed out. Note that the 
flow function and blockage factor were employed to 
calculate the new flow rate. Go to step 25. 

The original flow rate is used. The stator flow function 
and blockage factor based on the flow rate and measure- 
ment are printed. 

The display will read "INPUT ALPHA 1 (DEG), OR 0." 
Inputing any number other than 0 will fix the stator 
exit flow angle in degrees. This method assumes VU, 

1s calculated from the stator torque and adjusts the 
value of VA, to agree with the elected stator exit angle. 
By pressing "0 EXECUTE" the stator exit angle 1s not 
chosen. Go to step 27. 

By choosing to fix the stator exit flow angle, the 
value you selected will be printed. The value of K, 
which is the ratio of VA, determined by the momentum 
equation to VA, determined by continuity, will be 
printed. VAlues from 1.0 to 1.2 are within acceptable 


limits. Go to step 29. 


ee 





Zs 


ZG 


Zo « 


BU, 


Bil. 


BZ. 


oa. 


34. 


The value of the stator exit as calculated from 
measurements will be used. The display will read "ENTER 
K." Values between 1.0 and 1.2 are within acceptable 
limits. Press "(K) EXECUTE." This value will be 
printed and the program will CHAIN with "TTR 3." 

TTR 3 will evaluate Control Volume B and will CHAIN 
WLERMAETR 4." 

The display will read "CALC THEOR. LOSSES ? YES = l, 

NO = 0." This program option permits you to calculate 
both the theoretical rotor and stator loss coefficients. 
Copies of REF. 16 and 18 will be needed. By pressing 

"]) EXECUTE," you have elected to calculate the losses. 
Press "0 EXECUTE" if you do not desire to calculate 

the losses. Go to step 37. 

The display will read "DEL ALPHA = (SOME NUMBER) ." 

Note the number and press "CONTINUE EXECUTE." 

The display will read "BETA 1 = (SOME NUMBER)." Note 
the number and press "CONTINUE EXECUTE." 

The display will read "INPUT ZETA PO, FIG. 15, 1174 
VA?" iusto Fig. 15 in Re=:. 16. Select the value of 
ZETA PO which corresponds to the values of DEL ALPHA, 
BETA 1; and ALPHA EXIT. press ("ZETA PO) EXECUTE." 

The display will read "RE ROTOR = (SOME NUMBER) ." 

Note this number and press "CONTINUE EXECUTE." 

The display will read "INPUT K-RE, FIG. 23, GA1074VA2." 
Turn to Fig. 23 in Ref. 18 and select the value of Kon 
that corresponds to the rotor blade Reynolds number 
given in step 33. Press "(K-RE) EXECUTE." 


ee 





Bo. 


36. 


38. 


Bo. 


40. 


41. 


42. 


43. 


The display will read "M2 (IS) = (SOME NUMBER) ." 

Note this number and press "CONTINUE EXECUTE." 

The display will read "INPUT KM, FIG. 18, 1174va,." 
Turn to Fig. 18 in Ref. 16, and select the value of KM 
that corresponds to the blade exit Mach Number as 
calculated from isentropic conditions. Press " (KM) 
EXECUrE.” 

The program will CHAIN with "TTR 6." The display will 
read “ENTER NEXT RED REC. #." Input the same record 


number from step 12 and press "EXECUTE." 


The display will read "STORE RED. DATA ? YES i 

NO = 0." Press "1 EXECUTE" if you desire reduced data 
to be stored and used for plotting and tabulation. 
Press "Q EXECUTE" if you do not desire to store reduced 
data. Go to step 4l. 

You have elected to store the reduced data in the Mass 
Memory, and the words "REDUCED DATA IS STORED IN FILE/ 
REDDAT/RECORD (NUMBER FROM STEP 12)." will be printed. 
Electing not to store the reduced data the words "THIS 
DATA WAS NOT STORED" will be printed. 

The display will read "(CONT) (EXEC) MORE DATA - OR 
(END)ie 

If more data points are to be reduced, press "CONTINUE 
EXECUTE." The display will read "TAB RAW DATA ? YES = 
NO = 0." If you desire not to tabulate raw data, press 


"0 EXECUTE." Then the display will read "TAB REDUCED 


2 


i 





44. 


45. 


46. 


DATA ? YES = 1, NO = 0." If you desire not to tabulate 
reduced data, press "0 EXECUTE." By inputting zeros 

for both of these tabulation questions the program will 
display "REDUCE NEXT POINT," and will CHAIN with "TTR 

1."_ Go to step 12. Note: Tabulation is usually 
elected after all points are reduced for a run. 

To tabulate the raw data press "l EXECUTE." The program 
will GET "TTR 7." The display will read "EXECUTE RECORD 
#'S: LOWEST, HIGHEST." Input your first and last 
record numbers and press "EXECUTE." The display will 
read "OMIT 4 RECORDS. ENTER #'S OR 0'S." This option 
permits you to omit any points. Input the record 
numbers and press "EXECUTE." 

The display will read, read "REQUEST DATA CH #: 0, l 
OR 10." Press "1 EXECUTE" will cause the 48 measurements 
to be tabulated under the heading: TTR INPUT DATA. 
Afterwards, the calculator will request the next channel 
number. Press "10 EXECUTE" will cause all other channels 
to be tabulated and printed with the same heading. Next, 
press "0 EXECUTE." The display will read "TAB REDUCED 
DATA ? YES = 1, NO = 0." Press "1 EXECUTE" and the 
program will GET "TTR 8." Go to step 46. Press "0 
EXECUTE" and the program will GET "TTR 1." Go to step 
Ia 

By electing to tabulate the reduced data, the program 


will GET "TTR 8." The display will read "ENTER LOWEST, 


eS 





47. 


48. 


49. 


SUE 


HIGHEST REC. #." Input the first and last record 
numbers, and press "EXECUTE." 

The display will read "ENTER 1ST AND LAST PT. #." 
Input the two points, and press "EXECUTE. ' 

The display will read "ENTER RUN #." Input the run 
number, and press "EXECUTE." 

The velocities, Mach Number, and losses will be tabu- 
lated and printed. The tabulation will stop to allow 
the page to be removed from the printer. Then press 
-CONTINUOE EXECUTE. * 

Tabulation and printing will resume until "DATA TABULATION 
COMPLETE" is printed. The display will read "CHAIN 
WITH TAPE READER." There is a 3 second pause, and the 


program will GET "TTR 1." Go to step 12. 
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AO 
Al 
A3 
A4 
A5 
A8 
A9 
BO 
Bl 
BZ 
B3 
B4 
EBS 
B6 
B7 
B8 


BY 


C0 
Gal 


eZ 


TABLE B-2 


DEFINITION OF VARIABLES AND EQUATIONS 


Array 

Area of Control Volume (Fig. B-2, B-3) 

Area of Control Volume 

Alpha l 

Alpha 2 

Decision Variables 

Rotor Shroud Pressure Integrated Over Area 
Rotor Shroud Pressures Integrated Over Area 
FQ A (8) Ref. 5 

Beeeg mA. (16) Ret. 5 

Beta l 

Beta 2 

Units (2*CP*g) 

Rotor Shroud Pressures Integrated Over Area 
Rotor Shroud Pressures Integrated Over Area 
Rotor Shroud Pressures Integrated Over Area 
Del Alpha 


Mixing Loss Coefficient, Referred to Average Kinetic 
Energy, Eq. 7l, Ref. 16 


Array 
Gas constant (ft-lbf/1lbm-R) 
Specific Heat (BTU/1bm-R) 


Nozzle Diameter (in.) 


9 





eS 
C4 
eS 


C6 


C7 
C8 


C9 


DO 
Di 
D2 
D3 
D4 
D5 
D6 
D7 
D8 
D9 
E0 
BL 
E2 
E3 
B4 
E5 
E6 


E/ 


Gravitational Constant (32.174 Loni Useseae 
Pipe Diameter (upstream) 


Ratio (D, /D_) 


OZ =p 
Conversion factor: 0.16384 
¥2.036 


Discharge Coefficient 

Thermal Expansion Coefficient 

Compressibility Coefficient 

Array 

See Table B-l 

See Table B-l 

See Table B-l 

See Table B-l 

See Table B-l 

See Table B-l 

See Table B-l 

See Table B-l 

See Table B-l 

See Table B-l 

Effective Total-to-Static Efficiency Eq. A (18) Ref. 
Next record Number 

Input Variables 

See TRble B-l 

Profile Loss Coefficient of Turbine At Design Point 
Profile Loss Coefficient (G,,) Eq. 46, Ref. 16 
Secondary Flow Loss (2) Eqieaos,) Ref. ~ 16 


Mixing Loss Coefficient (C4) Eq. 74, Ref. 16 
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E8 


E9 


FO 
PL 
F2 
F3 
F4 
F5 
F6 


B/ 


F9 


GO 
Gl 
G2 
G6 
G7 
G8 


G9 


HO 
fants 
H2 


H5 


Stator Loss Theoretical (Total Blading Losses) 
Eq. 77, Ref. 16 


Rotor Loss Theoretical (Total Blading Losses) 
BG 7 7p Ret. LG 


Overall Axial Force 

Partial Sum of Overall Axial Force 

Pol*el 

Pa tot 

Stator Axial Force 

Closure Plate Force 

Ratio: (Thickness to Spacing) BG.873, REE. 16 


Mixing Loss Coefficient, Referred to Average Kinetic 
Energy, Eq. 71, Ref. 16 


Correction Factor for Blade Thickness, Eq. 47, Ref. 16 


Array 


y- 1 
vi 
Input Variable 


poele 

a ck 

Decision Variable 

Decision Variable 

Secondary Flow Loss, Eq. 63, Ref. 16 
Thickness to Spacing, Eq. 73, Ref. 16 
Pressure at Hub 

See Table B-l 

Horsepower 


Referred Horsepower 


Correction Factor For Effects of Reynolds Number 
(iernme) 
Re 


Zk 





H6 


H8 


Hy 


IQ 


fet 


I2 


3} 


I4 


i> 


I6 


lee? 


I8 


Eg 


J0 
J2 
J4 
J5 
J6 


J7 


Correction Factor For Blade Thickness (K,) Ba. 47, 
Ref. 16 


Density Appendix B, Ref. 5 
Reynolds Number of Rotor, Table B-2, Ref. 5 
First Point (Used in Tabulation Program) 


Ratio: Non-Dimensional Pressures At Various Tap 
hOcaL~ons In Sitator (Fig. 5) = DO in Table B-1 


Ratio: Non-Dimensional Pressures at Various Tap 
Locations In Stator (Fig. 5) - Dl in Table B-l 


Ratio: Non-Dimensional Pressures At Various Tap 
Locations In Stator (Fig. 5) - D2 in Table B-1l 


Ratio: Non-Dimensional Pressures At Various Tap 
Locations In Stator (Fig. 5) - D3 in Table B-l 


Ratio: Non-Dimensional Pressures At Various Tap 
Locations In Stator (Fig. 5) - D4 in Table B-1l 


Ratio: Non-Dimensional Pressures At Various Tap 
Locations In Stator (Fig. 5) - D5 in Table B-l 


Ratio: Non-Dimensional Pressures At Various Tap 
Locations In Stator (Fig. 5) - D6 in Table B-1l 


Ratio: Non-Dimensional Pressures At Various Tap 
Locations In Stator (Fig. 5) - D7 in Table B-l 


Ratio: Non-Dimensional Pressures At Various Tap 
Locations In Stator (Fig. 5) - D8 in Table B-l 


Ratio: Non-Dimensional Pressures At Various Tap 
Locations In Stator (Fig. 5) - D9 in Table B-l 


Mach # At Station 9 In Stator 

Joules Constant (778.16 Ft-lbf/BTU) 

See Table B-l 

Input Variable 

Profile Loss Coefficient Eq. 46, Ref. 16 
Velocity at Station 9 in Stator 


Non-Dimensional Stator Pressure 


nee 





J8 


J9 


KO 
Kl 
K2 
K3 
K4 
K5 
K6 
K7 
K8 


K9 


LO 
Ll 


2 


L3 


L4 


L5 


L6 


L7 


L8 


Velocity at EX-TIP 

Velocity at Station 5 in Stator 

Fourth Degree Polynomial For K.E Factor (Fig. 4) 
Viscosity (Sutherland Formula) 

Kinetic Energy Factor (G. E. Report) 

Beta (G. E. Report) 


K x B x AL x P. (KxK2xk6xQ0) 


L 
Isentropic Head Coefficient Eq. a (22), Ref. 5 
Slope of K. E. Curve (Newtonian Method) 

Area of Labyrinth 

K. E. Factor (Convergence) 

Second Degree Polynomial KfR,) For Nozzle (Fig. 5) 
See Table B-l 

Ce unt ez 

x- a“ + (xu2 - U2) Eq. A (21), Ref. 5 

Rotor Loss Coefficient Eq. A (21), Ref. 5 


Mach # At Various Pressure Tap Locations in Stator 
(Page. 5) 


Mach # At Various Pressure Tap Locations in Stator 
(Rag. 25) 


Mach # At Various Pressure Tap Locations in Stator 
CEG) 


Mach # At Various Pressure Tap Locations in Stator 
Ossi cpg Ss} 


Mach # At Various Pressure Tap Locations in Stator 
(Fig~ 95) 


Mach # At Various Pressure Tap Locations in Stator 
(Erone5) 


Mach # At Various Pressure Tap Locations in Stator 
CEarGie)) 
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L9 


MO 
M1 
M2 
M3 
M4 
M5 
M6 
M7 
M8 
M9 
NO 
N1 
N2 
N3 
N4 
N5 
N6 
N7 
N8 
N9 
00 
Ol 


O2 


os 


Mach # At Various Pressure Tap Locations in Stator 
(Fig. 5) 


Referred Moment 

Mach # At Station l 

Axial Mach # At Station 1l 
Mach # At Station 2 

Axial Mach # At Station 2 
See Table B-l 

See Table B-l 

Referred Stator Moment 
Mach # At Exit - Tip 

Mach # At SS - Aft 

See Table B-l 

See Table B-l 

See Table B-1l 

See Table B-l 

See Table B-l 

See Table B-l 

See Table B-1l 

See lane BS 

See Table B-l 

See Table B-l 

Incidence Loss Coefficient 
See Table B-l 


Velocities at Various Pressure Tap Locations In 
Stateone(Fig. 5} 


Velocities at Various Pressure Tap Locations In 
SiacOmme(rig. 5} 
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O04 


O5 


O6 


07 


08 


O9 


PQ 


Pl 


P2 


P3) 


P4 
PS 
P6 
17 
PY 
Q0 
Ql 
Q2 
Q3 
Q4 
@5 


Re 


Velocities at Various Pressure Tap Locations In 
Stator (Fig. 5) 


Velocities at Various Pressure Tap Locations In 
Stator (Fig: 5) 


Velocities at Various Pressure Tap Locations In 
Stator (Fig. 5) 


Velocities at Various Pressure Tap Locations In 
Stator (Fig. 5) 


Velocities at Various Pressure Tap Locations In 
Stacom tr ids. 5) 


Velocities at Various Pressure Tap Locations In 
Stator (Fig. 5) 


Reference Pressure 


(299 2 oat HG) 


P-HOOD/P-LAB (RATIO) 


W 
g Vt *Pt 


P,/P 
aon 
Po 


P,/Pt, 


Eq. A (10) Ref. 5 


2 
re 
vom ae awe 
eS es 


5 (Alpha Input) 
(sin a) 


(K Input) 


Pressure Ratio (Pro/Po) 


See Table 


Ba. 


Average Pressure At Rotor Exit 


See Table 
See Table 
2 ae 
See Table 


See Table 


Beer 


Bat 


Bele 


Bore 


Decision Variable 


Array 


Reynolds Number (Nozzle) 


25 





RO 


R2 
RS 
R4 
R8 


Ro 


SO 
Sl 
S2 
S3 
S4 
So 
S6 
S7 
S8 


S9 


TO 
T2 
3 
ES 
T6 
L 7 


T8 


Referred RPM (N//6) 

Reynolds Number (Labyrinth) 

Reva Am(4), Ret. 5 

Theoretical Degree of Reaction Eq. A (23), Ref. 


Actual Degree of Reaction Eq. A (24), Ref. 5 


See Table B-l 


Ratio (14.696/P_,_) 


Array 


Temperature Subroutine Variable 


oa 


Temperature Value 


Temperature Value 


See 


See 


See 


See 


See 


See 


See 


Pressure at Tip 


Reference Temperature 


Teo 


Y@ Temperature Ratio 


tn 


T 
Tep 
if 


S 


Table 
Table 
Table 
Table 
Table 
Table 
Table 


(Upstream Temperature °R) 


(Nozzle Temperature °R) 


(Plenum Temperature °R) 


Bae 
Bat 
Baa 
Bae 
Bae 
Ee 


B= 


(Fig. B-4) 


Ty (Fig. B-4) 





T (hig, B=-4) 


21s 
Array 
Mixing Loss Coefficient, Eq. 74, Ref. 16 


(Dy + D,)/48 


0 
(Dy + D.)/48 

Dimensionless Velocity 

Dimensionless Velocity 

RactOrmey fF.) 5/, Ret. 16 

Loss Due to Tip Clearance Flow, Eq. 75, Ref. 16 
Tip Clearance Loss Coefficient, Referred To Isen- 
tropic Conditions Eq. 76, Ref. 16 

Prof Py 

Stator Flow Function 

Array 


Tangential Rotor Velocity (Vea 


Limiting Velocity (V,_) 


12 Average Tangential Velocity Eq. A (11), Ref. 5 


al 
a2 


2 


See 


aL 
Eq. A (10), Ref. 5 
Stator Blockage Factor 
Array 
Total Flow Rate (W = Wee Wi ab) 


Labyrinth Flow Rate Based on K. E. Factor (G. E. 
Report) 


Za 


W2 
W3 
w4 
W5 
W6 
W7 
w8 


W9 


X0 
XL 
X1 
2 
X3 
X4 
X5 
X6 
X7 
X8 


X9 


YO 


pl: 


2 


Labyrinth Convergence Flow Rate 
Flow Rate (ASME Equation) 

Ratio (WO/g) 

Referred Flow Rate 


W, (See Fig. B-5) 


i 


W. (See Fig. B-5) 


2 


* = 
W 1 Wy 


(See Fig. B-5) 


Cos (By - 62) 


ae 


Array 
Dimensionless Velocity Xo 


Xo =n he (ve, Rei. 5 (K Input) 


Xai (Alpha Input) 
Xx = xX + X 2 
u2 a 


1/A PdA 


Xx Eq. A (7), Ref. 5 


u2 


X Eq. A (15), Ref. 5 


a2 
2 
Eee (Alpha Input) 
vacck—bjet 5 

Denominator of Eq. 63, Ref. 16 


Array 


For I0 
For Il 


For I2 
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wo 


Y4 


yD 


Y6 


ey 


Y8 


rg 


ZO 
ZL 
Z2 
Z3 
Z6 
Z7 


28 





p. ol 
ae ae ay 
£o 

p. xt 
ee (s— ae 
© 





Array 


Cl Eq. A (8a), Ref. 


FOr 


For 


For 


EOL 


For 


For 


For 


5 


Stator Loss Coefficient 


G72 EG.) 46116) ,7ReL. 5 


Rotor Loss Coefficient Eq. A (21), Ref. 5 


Non-Dimensional Pressure At Exit - Hub 


ES 


L6 


ay, 


I8 


ng 


Loss Due Tip Clearance, Eq. 


Rotor Loss Coefficient 
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(Other Factors) 





TABLE B- 3 
SIMPLE VARIABLE 


PROGRAM: 


XXKEXXKKKEKEEEX_EKEECEREE 
SPDR XPD IKD DT KPKPEP PPS 
REXEL XK SKEET i D<iebee RP 
RXR KD DPD PTR 
XXKKMT KKK KKK KEKE 
DKK XPD. Peb< Pox 
XK KL Deo PePePePh 
XKKRKREK KKK RK EKER RRR 
XK RPS 
SXRD DDD ePcx<pePeh< 
XXX | bx Tt | 
| DX 


XS 


XS 


XX | PDX 
s 


xX 


'D 


2 
2 


x 


SUBSCRIPT 


_ ali falaly|s|e|7jels 


2 


CC fe fc > 


Ax XX 
P: 
> 
x 
- 
= 
a 
xX 
oe 


SCORE SHEET: RECORD OF VARIABLES USED. 
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FIG. B-1 DATA REDUCTION SCHEMATIC 


TTR 1 
RAW DATA STORAGE 
RAW DATA PRINTED 
CHANNEL ASSIGNMENT 


TEMPERATURE SUBROUTINE 


TTR 2 
CALCULATE MASS FLOW RATE 


EVALUATION OF CONTROL 
VOLUME A ( FIG. B-2) 


ie 


EVALUATION OF CONTROL 
VOLUME B AND TEMPERATURES | 


TTR 4 


CALCULATE STATOR EXIT 
VELOCITIES (FIG. B=-5) 


CALCULATE THEORETICAL 
LOSSES 
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FIG. B-1 DATA REDUCTION SCHEMATIC 


TTR 5 


1.) PRINT REDUCED DATA 


TTR 6 


1.) STORE REDUCED DATA 


TTR 7 


1.) TABULATION OF RAW DATA 


TIR 8 


1.) TABULATION OF REDUCED DATA 





nae 








AIR /MLET 


rae 


Sf 








FIGURE B-2 CONTROL VOLUME a 


STATOR 


eee 6 ee oo ee 





FIGURE B-3 CONTROL VCLUME b 
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FIGURE B-4 THERMODYNAMIC PROCESS OF FLUIL 
IN AN AXIAL TURBINE STAGE 
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